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1 Introduction 

One of the most fascinating and unsolved problems of the theoretical physics 
of our century is the cosmological constant. Einstein introduced his cosmological 
constant A c in an attempt to generalize his original field equations. The modified 
field equations are 

R^u - 1 9^R + Kg^v = SttGT^, (1) 

where A c is the cosmological constant, G is the gravitational constant and T M „ is 
the energy-momentum tensor. By redefining 



jrtot _ rp _ J±c_ 



one can regain the original form of the field equations 

Rfiv — \ 9^vR = SnGT^J = 8irG (T^ + T^J) , (2) 

at the prize of introducing a vacuum energy density and vacuum stress-energy 
tensor 

A-r, ~nA 



8ttG ' " 

Alternatively, Eq. (1) can be cast into the form 



Pa - i T tiv - -Ph-g^v ■ 
be cast into the form, 

Rnu — 2 g^R + A-esgfiu = o , 
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where we have included the contribution of the vacuum energy density in the form 
Tp V = —(p)g^v In this case A c can be considered as the bare cosmological constant 

A cff = 8nGp cS = A c + 8ir G(p) . 

Experimentally we know that the effective energy density of the universe p c g is of 
the order 10~ 47 GeV 4 . A crude estimate of the Zero Point Energy (ZPE) of some 
field of mass m with a cutoff at the Planck scale gives 

This gives a difference of about 118 orders [1]. The approach to quantization of 
general relativity based on the following set of equations 



£zpe - =; I 7r4\/M^ « W 71 GeV 4 . (3) 



2 K G ijkl K ij TT kl ~^(R- 2A C ) 



* \9ij\ = (4) 



and 

-2Vi7r ij '*[ffy] = 0, (5) 

where R is the three-scalar curvature, A c is the bare cosmological constant and k — 
8irG, is known as Wheeler-De Witt equation (WDW) [2]. Eqs. (4) and (5) describe 
the wave function of the universe. The WDW equation represents invariance under 
time reparametrization in an operatorial form, while Eq. (5) represents invariance 
under diffeomorphism. Gijki is the supermetric defined as 

G ijkl = ^( 9ik9jl+9 u 9jk - 9ij9k i)- 
Note that the WDW equation can be cast into the form 



2KG ijkl ir ij n kl -^R 



2k 



*[ff«] = -— Ac*[ff«], 



which formally looks like an eigenvalue equation. In this paper, we would like 
to use the Wheeler-De Witt (WDW) equation to estimate (p). In particular, 
we will compute the ZPE due to massive and massless gravitons propagating on 
the Schwarzschild background. This choice is dictated by considering that the 
Schwarzschild solution represents the only non-trivial static spherical symmetric 
solution of the Vacuum Einstein equations. Therefore, in this context the ZPE 
can be attributed only to quantum fluctuations. The used method will be a vari- 
ational approach applied on gaussian wave functional. The rest of the paper is 
structured as follows, in section 2, we show how to apply the variational approach 
to the Wheeler-De Witt equation and we give some of the basic rules to solve such 
an equation approximated to second order in metric perturbation, in section 3, we 
analyze the spin-2 operator or the operator acting on transverse traceless tensors, 
in section 4 we use the zeta function to regularize the divergences coming from the 
evaluation of the ZPE for TT tensors and we write the renormalization group equa- 
tion, in section 5 we use the same procedure of section 4, but for massive gravitons. 
We summarize and conclude in section 6. 
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2 The Wheeler— De Witt equation and the cosmological constant 



The WDW equation (4), written as an eigenvalue equation, can be cast into the 
form 

A s * [ 9lJ ] = -A'(aO* [9ij] , (6) 

where 

A E = 2 K G ijk i7r i ^ kl - y^R, 



2k 



K 



We, now multiply Eq. (6) by \P* [gij] and we functionally integrate over the three 
spatial metric g^, then after an integration over the hypersurface S, one can for- 
mally re-write the WDW equation as 



V 



1 

V 7 



is d3a;A s 



<*!*> 



= A'. 



(7) 



The formal eigenvalue equation is a simple manipulation of Eq. (4). However, we 
gain more information if we consider a separation of the spatial part of the metric 
into a background term, , and a perturbation, hij , 



gij gij ~l~ hij . 



Thus eq. (7) becomes 



il- 



ls d 3 * 



A< 3 0) +A^+A£' + .. 



;(2) 



<*!*> 



A'* [g 



(8) 



where A^ represents the i th order of perturbation in hij. By observing that the 
kinetic part of Ag is quadratic in the momenta, we only need to expand the three- 
scalar curvature J d^x^/gR^ 3 ' up to quadratic order and we get 



d 3 xVff 



i ftAft + ± h li Ahu - i h ij ViVM + 



+ i hViVih" - \ h ij R ia h^ + § hRijh ij + \h ft 



(9) 



where ft is the trace of hij and R^ is the three dimensional scalar curvature. To 
explicitly make calculations, we need an orthogonal decomposition for both mj and 
to disentangle gauge modes from physical deformations. We define the inner 
product 

(h,k) := / ^gG ijkl hij(x)k k i(x)d 3 x, 



by means of the inverse WDW metric Gijki, to have a metric on the space of 
deformations, i.e. a quadratic form on the tangent space at hij, with 

Qijkl = ( g ik g jl + g il g jk _ 2 gVg kl ) . 
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The inverse metric is defined on cotangent space and it assumes the form 
(P, V) ■= \ 7 I<-. • {x)q U (x) d 3 x , 

so that 

G^ nm G nmkl = i fa' + 5j5i) ■ 

Note that in this scheme the "inverse metric" is actually the WDW metric de- 
fined on phase space. The desired decomposition on the tangent space of 3-metric 
deformations [3-6] is: 

In, = Ihgij + iL^+h^, (10) 
where the operator L maps into symmetric tracefree tensors 

(i€)« = v^ + v^-ifl«(v-0- (ii) 

Thus the inner product between three-geometries becomes 

(h,h) := / ^gG i ^ kl h ij (x)h kl (x)d 3 x = f ^g[-^h 2 + (L0 ij (L^) ij +h ii± hj ] ]. 
JE je 

_ ^ (12) 

With the orthogonal decomposition in hand we can define the trial wave functional 
as 

* [hn ( <?)] = ^ H {*)} * [hi &)] * [hrr ( * )] - (13) 



where 



* [4 (1?)] = cxp {-i <(L£) (L0)!,J , 

vf (1?)] = cxp {-I (hK-'h)^ e } . 

The symbol "_L" denotes the transverse-traceless tensor (TT) (spin 2) of the pertur- 
bation, while the symbol "||" denotes the longitudinal part (spin 1) of the pertur- 
bation. Finally, the symbol "trace" denotes the scalar part of the perturbation. J\f 
is a normalization factor, (•, -) x . y denotes space integration and K^ 1 is the inverse 
"propagator''' . We will fix our attention to the TT tensor sector of the perturbation 
representing the graviton. Therefore, representation (13) reduces to 

* [hij ("?)] = AAcxp j-i (hK-'h^Jf . (14) 

Actually there is no reason to neglect longitudinal perturbations. However, follow- 
ing the analysis of Mazur and Mottola of Ref. [5] on the perturbation decomposition, 
we can discover that the relevant components can be restricted to the TT modes 
and to the trace modes. Moreover, for certain backgrounds TT tensors can be a 
source of instability as shown in Rcfs. [7]. Even the trace part can be regarded as 
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a source of instability. Indeed this is usually termed conformal instability. The ap- 
pearance of an instability on the TT modes is known as non conformal instability. 
This means that does not exist a gauge choice that can eliminate negative modes. 
To proceed with Eq. (8), we need to know the action of some basic operators on 
"J [hij]. The action of the operator h^ on \&) — * [hij] is realized by [8] 

M*)l*) = hij (2)* [hij]. 
The action of the operator mj on | , in general, is 

S 



I*) 



Shij (a?) 



* [hij] , 



while the inner product is defined by the functional integration: 
(* i | *2> = J Phij] *I [hij] * 2 [hki] ■ 



We demand that 



k d 3 zA E 



V 



(*!*) 



V 



fV[ gi j]** [hij]* [hi 



(15) 



be stationary against arbitrary variations of * [hij]. Note that Eq. (15) can be 
considered as the variational analog of a Sturm-Liouville problem with the cosmo- 
logical constant regarded as the associated eigenvalue. Therefore the solution of 

Eq. (7) corresponds to the minimum of Eq. (15). The form of (ty A E vf^ can be 

computed with the help of the wave functional (14) and with the help of 

( (* |M"?)|*) 



(*\hjjjx)hki {It)]*) 
(*|*) 



= KijM (a?,"y ) 



Extracting the TT tensor contribution, we get 



ijkl 



2kK- 1± (x 1 x)ijki + — (A 2 )° K^(x, x) 



2k 



iakl 



The propagator K ± (x, x) iaki can be represented as 



2A(r) 



(16) 



(17) 



where h\^ >± (a? ) are the eigenfunctions of A 2 . r denotes a complete set of indices 
and A(t) are a set of variational parameters to be determined by the minimization 
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of Eq. (16). The expectation value of is easily obtained by inserting the form 
of the propagator into Eq. (16) 



4 

T 2 — 1 



2«A,(r) ■ ^ (T) 



2/cAi(r) 



By minimizing with respect to the variational function Aj(r), we obtain the total 
one loop energy density for TT tensors 



4 

r 

The above expression makes sense only for tof(r) > 0. 



(18) 



3 The transverse traceless (TT) spin 2 operator for the Schwarzschild 
metric and the W.K.B. approximation 

The spin-two operator for the Schwarzschild metric is defined by 

{A 2 h^)\ := - (A T /> TT )J + 2 (Rh^)\ , (19) 

where the transverse-traceless (TT) tensor for the quantum fluctuation is obtained 
by the following decomposition 

This implies that (h T )'' i <5* = 0. The transversality condition is applied on (h T y. 
and becomes Vj (^ T )^ = 0. Thus 



_(A T ^)J = -A fl (^)J + « (l-^) 



(20) 



where Ag is the scalar curved Laplacian, whose form is 



A -(l 2MG \ d2 | f 1r-ZMG \ d L 2 
S \ r ) dr 2 \ r 2 ) dr r 2 

and Rj is the mixed Ricci tensor whose components are: 

2MG MG MG\ 

Ri — { — 



I' 

This implies that the scalar curvature is traceless. We are therefore led to study 
the following eigenvalue equation 

(A^Wfcj, (22) 
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where ui 2 is the eigenvalue of the corresponding equation. In doing so, we follow 
Regge and Wheeler in analyzing the equation as modes of definite frequency, an- 
gular momentum and parity [9]. In particular, our choice for the three-dimensional 
gravitational perturbation is represented by its even-parity form 



(h e ™)) (r, <?, 0) = diag [H(r),K(r), L(r)} Y lm {#, 0) . 



(23) 



with 



H{r) = h\{r)-\h{r), 

K{r) = hl{r)-\h{r), 

L{r) = hl{r)-\h{r). 

From the transversality condition we obtain h\{r) = h\{r). Then K(r) = L(r). For 
a generic value of the angular momentum L, representation (23) joined to Eq. (20) 
lead to the following system of PDE's 



2MG\ AMG 



-A 



S "I" -o 



2MG\ 2MG\ 



Defining reduced fields 



H(r) = 



Mr) 



K(r) = 



H(r) = ^ u H(r) 
j K(r) = ^ u K{r) . 

Mr) 



(24) 



and passing to the proper geodesic distance from the throat of the bridge 

dr 



dx = ± 



1 



2MG 



(25) 



the system (24) becomes 



d 2 



h{x) =ul l f 1 (x), 

f 2 (x) = uhM x ) 



(26) 



with 



Vl{r) = l l±ll + Ulir) . 
V 2 (r) = l -V±^ + U 2 (r). 



where we have defined r = r(x) and 



\Ui{r) = 


' 6 




3MG" 
















{ U 2 (r) = 


' 6 




3MG" 
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Note that 



f/i(r)>0, when r > § MG , 

C/i(r) < when 2MG < r < | MG , 



(27) 



C/ 2 (r) > 



Vr e [2MG,+oo) . 



The functions U\{r) and ?/2(^) play the role of two r-dependent effective masses 
m\(r) and m\(r), respectively. In order to use the WKB approximation, we define 
two r-dependent radial wave numbers k\ (r,l,ui in i) and k 2 {r,l,u)2, n l) 



k\ (r, I, u hn i) = u? l nl - X) - m\(r) , 

k\ {r,l,uj 2 ,nl) = u\ nl - _ m 2( r ) 



for r > f MG. When 2MG < r < § MG, fc? (r, l,<J hnl ) becomes 



k\ (r,l,(Ji, n i) = Lo\ nl 



(28) 



(29) 



4 One loop energy regularization and renormalization 

In this section, we proceed to evaluate Eq. (18). The method is equivalent to 
the scattering phase shift method and to the same method used to compute the 
entropy in the brick wall model. We begin by counting the number of modes with 
frequency less than u>i, i = 1,2. This is given approximately by 



ff(w i ) = 5^i/ i (i,w i )(2i + l), 



(30) 



where Vi (I, cji), i = 1, 2 is the number of nodes in the mode with (l,u>i), such that 
(r = r(x)) 

1 f +co / 

Vi(l,Wi) = —l dxykf(r,l,uji). (31) 



Here it is understood that the integration with respect to x and / is taken over 
those values which satisfy k\ (r, > 0, i = 1,2. With the help of Eqs. (30, 31), 
we obtain the one loop total energy for TT tensors 



2 A z' + oc r r+oo 

to&J-oo Jo 



duji 



du>i 



By extracting the energy density contributing to the cosmological constant, we get 

(32) 



A = Ai + A 2 = p\ + p 2 = 



K 



16tt 2 



o 



lj\ \Ju)\ — m\{r) du)i 



16tt 2 



o>2 — m|(r) dw 2 , 
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where we have included an additional 4ir coming from the angular integration. We 
use the zeta function regularization method to compute the energy densities pi and 
p 2 - Note that this procedure is completely equivalent to the subtraction procedure 
of the Casimir energy computation where the zero point energy (ZPE) in different 
backgrounds with the same asymptotic properties is involved. To this purpose, we 
introduce the additional mass parameter p in order to restore the correct dimension 
for the regularized quantities. Such an arbitrary mass scale emerges unavoidably 
in any regularization scheme. Then we have 



Pi(e) 



1 



,2e 



16tt 2 



+00 



dujj 



(w?-m?(r)) 



e-l/2 ' 



(33) 



The integration has to be meant in the range where ujf — mf(r) > 1 ). One gets 



I + i„ 

e 



P 



+ 21n2- - 
2 



(34) 



2567T 2 [e \ m f( r ) 
i = 1, 2. In order to renormalize the divergent ZPE, we write 
A = 8nG( Pl (e) + p 2 (e) + pi (p) + p 2 (p)) , 

where we have separated the divergent part from the finite part. To handle with 
the divergent energy density we extract the divergent part of A, in the limit e — > 
and we set 



A 



div 



G 
32tt£ 



(mi(r)+m 4 2 (r)) 



Thus, the renormalization is performed via the absorption of the divergent part 
into the re-definition of the bare classical constant A 

A -> A + A div . 

The remaining finite value for the cosmological constant reads 



1 



Ao 

8ttG 256tt 2 



m\(r) 



In 



P 



\m\(r)\ 



+ 2 In 2 - - 
2 



+ 



In 



m\(r) 



2 In 2 - - 
2 



(35) 



(Pi(^) + P2(p)) = pjs (M, r) . 



The quantity in Eq. (35) depends on the arbitrary mass scale p. It is appropriate 
to use the renormalization group equation to eliminate such a dependence. To this 
aim, we impose that [10] 



1 dAfi T (p) 
P ' 



8ttG 



dp 



d TT / \ 

^dp PcS ^' r ' ' 



(36) 



1 ) Details of the calculation can be found in the Appendix. 
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Solving it we find that the renormalized constant Ao should be treated as a running 
one in the sense that it varies provided that the scale m is changing 



Ao(m,»") = A (Mo,r) 



G 

167T 



(m*(r)+m*(r))ln 



Mo 



(37) 



Substituting Eq. (37) into Eq. (35) we find 



Ao(^ ,r) 
8ttG 



1 



256tt 2 
+mi(r) 



m\{r) 



In 



In 



Mo 



\m\ (r)| 
Mo 



21n2 



21n2+ - 



+ 



(38) 



In order to fix the dependence of A on r and M, we find the minimum of Ao(mo, r )- 
To this aim, last equation can be cast into the form 2 ) 



Aq(mq,7') 
8ttG 



Mo 



256tt 2 



x 2 {r) 



In 



\x(r)\ 



1 

+ 2 



y 2 (r) 



In 



(39) 



where x(r) = ±m 2 (r)/ij, 2 l and y(r) = m|(r)/^g. Now we find the extrema of 
A (mo! x ( r ), y( r )) m the range § MG < r and we get 



x(r) = 0, 
y(r)=0, 



which is never satisfied and 



x(r) = 4/c , 
y(r) = 4/c , 



m i( r ) = 4Mo/ e : 
m\{f) = 4//g/e, 



(40) 



(41) 



which implies M = and r — \/3e/2Mo- On the other hand, in the range 2MG < 
r < | MG, we get again 

'ir(r) =0, 
y(r) = 0, 



which has no solution and 



which implies 



-m\{r) = 4^/e, 
. ml{r) = 4^o7e, 

M = 4 / u r 3 /3eG, , 
f = V6c/4mo • 



(42) 



(43) 



2 ) Recall Eqs. (28,29), showing a change of sign in m\ (r). Even if this is not the most appropriate 
notation to indicate a change of sign in a quantity looking like a "square effective mass", this 
reveals useful in the zeta function regularization and in the search for extrema. 
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Eq. (39) evaluated on the minimum, now becomes 

Ao(M,r)=fg,. (44) 

It is interesting to note that thanks to the renormalization group equation (36), 
we can directly compute Ao at the scale and only with the help of Eq. (37), we 
have access at the scale /x. 

5 One loop energy regularization and renormalization for massive 

gravitons 

The question of massive gravitons is quite delicate. A tentative to introduce a 
mass in the general framework has been done by Boulware and Deser [11], with 
the conclusion that the theory is unstable and produces ghosts. However, at the 
linearized level the Pauli-Fierz term [12] 



>P.F. 



m 2 f 

I & 4 x^—g [Wh^ - h 2 ] , (45) 



does not introduce ghosts. m g is the graviton mass. Following Rubakov [13], the 
Pauli Fierz term can be rewritten in such a way to explicitly violate Lorentz sym- 
metry, but to preserve the three-dimensional Euclidean symmetry In Minkowski 
space it takes the form 

Sm = —^-[ d 4 x^[m 2 h 00 h Q0 + 2m 2 h 0t h Ql -m 2 h^h l] + 

ok J M (46) 

• ///;;/'"//,, - 2m 2 h 00 h u ] . 

A comparison between the massive action (46) and the Pauli-Fierz term shows that 
they can be set equal if we make the following choice 3 ) 

m 2 , = , m\ = m 2 = to 2 = m 2 = m 2 > . 
If we fix the attention on the case 

m§ = m\ = to 2 = m 2 = , m 2 . = m 2 > , (47) 
we can see that the trace part disappears and we get 

w 2 f , i r ■ n 

Sm = ~d J dx ^i hVh v]- 

The corresponding term in the linearized hamiltonian will be simply 



TTl 2 f 

H m = —^J d 3 xN^ [Who] . 



3 ) See also Dubovski [14] for a detailed discussion about the different choices of mi, m,2, m-3 
and m4 
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This means that Eq. (19), will be modified into 



(A 2 h^y t := - (A T h^ + 2 (Rh^)\ + K^ TT ) 



■ TT\J 



TT\3 



2. TT\J 



(48) 



Therefore, the square effective mass will be modified by adding the term m 2 g . Note 
that, while m?,(r) is constant in sign, m 2 (r) is not. Indeed, for the critical value 



5MG/2, m\{f) 



m. 



and in the range (2MG, 5MG/2) for some values of 



nig, m\(f) can be negative. It is interesting therefore concentrate in this range. 
To further proceed, we observe that m\(r) and m|(r) can be recast into a more 
suggestive and useful form, namely 

j m\(r) = m 2 g + U\{r) = m 2 g + m\{r,M) - m|(r,M) , 
\m 2 2 (r) =m 2 g + U 2 {r) = m 2 g + mf(r, M) + m|(r, M) , 

where mf(r,M) -» when r -» oo or r -» 2MG and m|(r,M) = 3MG/r 3 . 
Nevertheless, in the above mentioned range m\(r, M) is negligible when compared 
with m 2 (r, M). So, in a first approximation we can write 



m\ (r 



](r ,M) =m 2 -m 2 (M), 



m\(r) ~ TOg + mKro, M) — m 2 + m§(M) . 



where we have defined a parameter ro > 2MG and mj)(M) = 3AfG/rQ. The main 
reason for introducing a new parameter resides in the fluctuation of the horizon 
that forbids any kind of approach. Of course the quantum fluctuation must obey 
the uncertainty relations. Thus, the analogue of Eq. (38) for massive gravitons 
becomes 



8ttG 



256^ U m » 



i 2 (M)y 



Til, 



In 



m 2 (M)Y 



In 



m 2 g - ml(M)\ 
Mo 



21n2+ - 



m 2 +ml{M) 



Mo 



21n2- 



+ 



(49) 



We can now discuss three cases: 



1. m 2 > ml{M), 



2. m 2 = m 2 {M), 



3. m 2 < ml(M). 
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In case 1), we can rearrange Eq. (49) to obtain 



A (^ ,r) 
8ttG 



1 



2(m4 + TO 4 (M) ) 



256tt 2 
+ +(m^ + mJ(M)) 



In'™* 



In 



+2m*m z (M) 

,4 



In 1 + 



l-#) + ]nf 
to 2 J V 



mg(M) 



+ 
1- 

2/ 



T7Jq ( Af 



m„ 



256tt 2 



2 In 



The last term can be rearranged to give 



128vr 2 



In I ^ 



Art 



M 



where we have introduced an intermediate scale defined by 



At = Mo ex P 



3mg(M) \ 
2m4 



(50) 



With the help of Eq. (50), the computation of the minimum of Aq T is more simple. 
Indeed, if 

.2 



m„ 



x = 



A becomes 



A . M (iio,x) = - 



In a; + 



(51) 



As a function of x, A ,M(tto,x) vanishes for x = and x — e x l 2 and when 
x G [0, cxp (—5)] , Aj"^- (/io, a;) > 0. It has a maximum for 



1 

x = - 
e 



2 4^4 4^ 2 



cxp - 



3mg(M) 
2m4 



and its value is 



or 



Ao,m(Mo,z) = n __ = ex P 



27re 2 2?re 2 
G 

A ,m(mo,s) = ^: m ff ex P 



3m^(M) 



Hi, 



3m^(M) 
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In case 2), Eq. (49) becomes 



AoOo,r) A (^o) 



8ttG 



8ttG 



128tt 2 




or 



Ao( M o) _ mg(M) 
SttG 128tt 2 

Again we define a dimensionless variable 



In 



m„ 



4 M 2 J + 2 



x = 



and we get 



A 0! o(^o,a;) 
8ttG 



lnx+- 



(52) 



The formal expression of Eq. (52) is very close to Eq. (51) and indeed the extrema 
are in the same position of the scale variable x, even if the meaning of the scale 
is here different. Ao,o(/-to,x) vanishes for x = and x = 4e -1 / 2 . In this range, 
^0.0(^0 7 x) > and it has a minimum located in 



x = 



m l = 



(53) 



and 



or 



A ,o(Mo,x) = 



G 



A ,o(Mo,x) = T^rrig 



2ire 2 
G 

~ 32tt 



Finally the case 3 ) leads to 
A (Mo,r) m^M) 



8irG 



256tt 2 



2 In 



) + + 



m 2 (M) 



The last term can be rearranged to give 

mj(M) f ( ml{M) \ I 
128^ 2 [ \ ) + 2 

where we have introduced another intermediate scale 



2 2 ( 3m s 
M ro = Mo ex P -7 



2m^(M) 
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By repeating the same procedure of previous cases, we define 

_ ml(M) 



and we get 

Ao, m (Mo,a;) 

Also this case has a maximum for 



G/x 4 



m x 2 



\nx+- 



1 

x = - 

e 



mn(M) = = cxp 

e e 



3m„ 



2m 4 (M) 



and 



or 



3to 4 



(54) 



Ao,mM = = gjexp ^-pjy y 

G 4 / 3m 4 „ \ 

A ,m(M0,S) = — mo (M)exp^^J. 

Remark. Note that in any case, the maximum of A corresponds to the minimum 
of the energy density. 



6 Summary and conclusions 

In this paper, we have considered how to extract information on the cosmolog- 
ical constant using the Wheeler Dc Witt equation when the graviton is massless 
and massive. In particular, by means of a variational approach and a orthogonal 
decomposition of the modes, we have studied the contribution of the transverse- 
traceless tensors in a Schwarzschild background. The use of the zeta function and 
a renormalization group equation have led to three different cases: 

' to 2 » m 2 (M) f A ,m(Mo, x) = G^J (2irc 2 ) exp (-3m 4 (M)/m 4 ) 

< to 2 = m 2 Q (M) =► J A ,o(Mo,S) = Gf4/ (2tc 2 ) 
k ml « m 2 (M) { Ao, m ( Mo , x) = G^/ (27re 2 ) cxp (-3m 4 /m 4 (M)) 

As we can see, the case "extreme" , where the graviton mass is completely screened 
by the curvature "mass" seems to have the biggest value. We recall that the high- 
est is the value of Ao , the lowest is the value of the energy density. However, 
the expression of the extreme case coincides with the mass-less graviton discussed 
in section 4. In that paper, it is the curvature "mass" which plays the role of 
the mass of the graviton and contributes to the cosmological constant. So it ap- 
pears that the gravitational field in the background of the Schwarzschild metric 
generates a "mass" term, because of the curvature and this term disappears when 
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the Schwarzschild mass goes to zero. This leads to the conclusion that fluctuations 
around Minkowski space do not create a cosmological constant in absence of matter 
fields. Nevertheless, this behavior works if we accept that near the throat, vacuum 
fluctuations come into play forbidding to reach the throat itself. If this is not the 
case and the throat can be reached, then the curvature "moss" becomes completely 
non-perturbative when the Schwarzschild mass M — > 0. If we choose to fix the 
renormalization point /z = m p , we obtain approximately Aq"(M, f) ~ 10 37 GeV 2 
which, in terms of energy density is in agreement with the estimate of Eq. (3). Once 
fixed the scale ^o, we can see what happens at the cosmological constant at the scale 
/i, by means of Eq. (37). What we see is that the cosmological constant is vanishing 
at the sub-planckian scale fi = m p e -1 / 4 , but unfortunately is a scale which is very 
far from the nowadays observations. Note that, because of the condition (53), the 
graviton mass becomes proportional to the "Planck mass" , which is of the order 
10 16 GeV, while the upper bound in cV is of the order Kr 24 -l(r 29 eV [15]. A 
quite curious thing comes on the estimate on the "square graviton mass" , which in 
this context is closely related to the cosmological constant. Indeed, from Eq. (53) 
applied on the square mass, we get 

m 2 oc M 2 ~ 10 32 GeV 2 = 10 50 cV 2 , 

while the experimental upper bound is of the order 

(ml) oc 1(T 48 - 1CT 58 cV 2 , 

which gives a difference of about 10 98 - 10 108 orders. This discrepancy strongly re- 
call the difference of the cosmological constant estimated at the Planck scale with 
that measured in the space where we live. However, the analysis is not complete. 
Indeed, we have studied the spectrum in a W.K.B. approximation with the follow- 
ing condition fc?(r, Z,Wj) > 0, i = 1, 2. Thus to complete the analysis, we need 
to consider the possible existence of nonconformal unstable modes, like the ones 
discovered in Refs. [7]. If such an instability appears, this does not mean that we 
have to reject the solution. In fact in Ref. [16], we have shown how to cure such 
a problem. In that context, a model of "space-time foam" has been introduced in 
a large N wormhole approach reproducing a correct decreasing of the cosmological 
constant and simultaneously a stabilization of the system under examination. Un- 
fortunately in that approach a renormalization scheme was missing and a W.K.B. 
approximation on the wave function has been used to recover a Schr6dinger-like 
equation. The possible next step is to repeat the scheme we have adopted here in 
a large N context, to recover the correct vanishing behavior of the cosmological 
constant. 
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A The zeta function regularization 



In this appendix, we report details on computation leading to expression 
We begin with the following integral 



P(e) = { 



1+ = /j; 



+00 



duj- 



(uj 2 + m 2 (r)) e - 1/2 ' 

+oc 2 £-1/2 

2e ' dec- 



/" 

J 



10 



(uj 2 — m 2 (r)) 



with m 2 (r) > 0. 



A.l J_|_ computation 

If we define t — lo/ yjm 2 {r), the integral 1+ in Eq. (55) becomes 



"+00 



p{s) — fi m (r) 



-_^_^-*(,)B (§,.-, 



m 2 (r)J r(e-i)' 
where we have used the following identities involving the beta function 

£2x-l 

Re x > , Re y > 



/•+OO 

B(x,y) = 2 dt- 
Jo ( 



( t 2 + if+v 

related to the gamma function by means of 

B(x y) T{x)r(y) 
Taking into account the following relations for the T-function 

r(i + £ ) 



r( £ -2) = 
r( £ -i) = 



£ ( £ -l)( £ -2)' 

r(e + |) 



and the expansion for small £ 

r(l+e) = l-7£ + 0(e 2 ), 
r(e + i) = r(i) - er(i)( 7 + 2 In 2) + O (e 2 ) , 
x £ = l + eln.T + (9(e 2 ) , 

where 7 is the Euler's constant, we find 



Pie) = 



m 4 (r) 
16 



1 



+ ln 



n 2 \ 1 
+ 21n2- - 
m 2 r / 2 
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A. 2 I- computation 

If we define t = to/ \Jm 2 (r), the integral J_ in Eq. (55) becomes 
p(s) = ^m^(r) £°° dt ^_^ e _ 1/2 = 1 ^m 4 - 2£ (r)B ( e - 2, | - s 



1 M^(r) r ^-^- 2 > = - * mV) fd^V r - e ) r ( e - 2), 



2^ vv r(-|) " 4^ v ^m 2 (r) 

where we have used the following identity involving the beta function 

1 b(i-v-£,J) = [ + °° dtt*- 1 ^-!)"- 1 



P \ P 

p > , Re ^ > , Reri < p — p Re f 

and the reflection formula 

T(z)T(l - z) = -zT(-z)T(z) . 
From the first of Eqs. (56) and from the expansion for small e 

T(| -e) = r(|)(l - e(- 7 - 21n2 + 2)) + O (e 2 ) 



x e = l + elnx + O (e 



2\ 



we find 



-+ln + 21n2-- 

16 \_s \m z (r) ) 2 
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